o 



3 



Stochastic particle annihilation: a model of state reduction in 
relativistic quantum field theory 

Daniel J. BedinghairQ 

Blackett Laboratory, Imperial College, London SW7 2BZ, UK 
(Dated: February 9, 2008) 

A model of state reduction in relativistic quantum field theory involving a nonlin- 
ear stochastic extension of Schrodinger's equation is outlined. The eigenstates of the 
annihilation operator are chosen as the preferred basis onto which reduction occurs. 



These are the coherent states which saturate the bound of the Heisenberg uncer- 



tainty relation, exhibiting classical-like behavior. The quantum harmonic oscillator 
is studied in detail before generalizing to relativistic scalar quantum field theory. 
The infinite rates of increase in energy density which have plagued recent relativistic 
proposals of dynamical state reduction are absent in this model. This is because 
the state evolution equation does not drive particle creation from the vacuum. The 
£NJ ■ model requires the specification of a preferred sequence of space-like hyper-surfaces 

supporting the time-like state evolution. However, it is shown that the choice of 
preferred surfaces has no effect on perturbative results to second order in the cou- 
pling parameter. It is demonstrated how state reduction to a charge density basis 
can be induced in fermionic matter via an appropriate coupling to a bosonic field 
' undergoing this mechanism. 
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~$ I. INTRODUCTION 

Much of the peculiar behavior associated with quantum physics results from the fact 
that, although a quantum system can be in a superposition of different states, whenever we 
make measurements involving macroscopic apparatus, a definite state is always registered. 
The transition from a superposition to a definite state is not described by Schrodinger's 
equation. How then, if the constituents of the apparatus are also described by Schrodinger's 
equation, does this quantum state reduction come about? 

Stochastic generalizations of Schrodinger's equation have been proposed by a number of 
authors in answer to the problem of measurement [j], 0, 0, 0, 0] (for a review see 0, 0] ) . The 
key idea is that measurement is understood as the realization of a random process in the 
Hilbert space of state vectors where unwanted superpositions of states are unstable. The 
appeal of these models rests on two fundamental properties: (i) they reproduce quantum 
effects on small scales with negligible modification to standard quantum theory, and (ii) 
they lead to rapid, objective state vector collapse on large scales with probabilities given 
by the laws of standard quantum mechanics. The result is that superpositions of states for 
macroscopic objects are suppressed whilst individual particles continue to behave according 



*Electronic address: d.bedingham@imperial.ac.uk 



2 



to quantum theory. 

The usual approach is to substitute Schrodinger's equation with a quantum state diffusion 
equation of the form 

d\<p t ) = (Cdt + A-dX t )\<f )t ). (1) 

Here {X t } is a (vector-valued) Ito process and A, C are operators (the Schrodinger equation 
can be recovered by setting C = —iH and A = 0). With appropriate choices for the drift 
and volatility of {X t } the quantum state typically evolves into an eigenstate of the operator 
A. The choice of A leads to a preferred basis. In the quantum mechanical case, the standard 
choice is a locally averaged position state basis in order to reproduce the definite localization 
of objects at the classical scale. Another idea is to use an energy state basis @, @, These 
models have the desirable property that energy is conserved in expectation. A general 
solution to the energy-based state diffusion with time-dependent coupling has recently been 
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At present, non-relativistic proposals are seen to have sufficiently negligible effects on 
the quantum scale in order to be indistinguishable from standard quantum theory for cur- 
rent experimental technologies [13]. At the same time these proposals offer a consistent 
understanding of classical and quantum domains. However, so far, relativistic field theoretic 
formulations generally predict an infinite rate of particle creation due to the coupling of a 
classical stochastic field to a quantum scalar field @, [H, 14, 15]. Some previous attempts 



to resolve this problem have involved modifying the stochastic field to prevent high-energy 
excitations 0, EE], or coupling the noise source not locally to the quantum field but to the 
integral of quantum fields over some space-time region [l7f. A quantum mechanical model 
for a relativistic particle has been developed in reference [18| although this model does not 
include interactions. 

In this paper we outline an alternative proposal in which the stochastic field is coupled 
only to the annihilation operators of the quantum scalar field (via a local interaction term). 
The scalar field cannot then be excited by the stochastic field. As a consequence the infinite 
rates of energy increase are avoided. Instead we see an expected energy loss to the stochastic 
field which can be controlled to a negligibly small level by an appropriate choice for the 
coupling parameter. A related idea has been employed in reference [19| to control energy 
increase in models of non-relativistic state reduction. 

We find that in order to construct a satisfactory model of state reduction in relativistic 
quantum field theory we must assume a preferred sequence of space-like hyper-surfaces 
supporting the evolution of the quantum state. The reason is that the stochastic field 
is coupled to local operators which do not commute at space-like separation. The state 
evolution equation is therefore path-dependent. The fixed sequence of space- like hyper- 
surfaces constrains the evolution such that only one path is possible, ensuring a well defined 
evolving state. We do not propose a rule for how the surfaces are chosen and regard them 
hidden property of the state. 

Our state evolution equations are of relativistically invariant form so that all observers 
will agree on outcomes. However, the choice of surface is responsible for identifying a 
preferentially selected local frame. The idea that dynamical reduction models might break 
Lorentz invariance in this way has been suggested before by Pearle [2(j, who considered 
a stochastic field coupled to a generalized mass-density field which does not commute at 
space-like separation. There it was shown that the commutator decays on a length scale 
corresponding to the particle's Compton wavelength, providing a sense in which the model 
is quasi-relativistic. 
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By performing perturbative calculations involving an expansion in the coupling parameter 
we are able to quantify the effect of a particular choice of the space-like hyper- surfaces. We 
find that the choice has no effect on the lowest order expressions describing state reduction. 
This offers an alternative way to understand the quasi-relativistic nature of this type of 
model. 

We will see that the quantum state evolves towards the eigenstates of the annihilation 
operators. In quantum mechanics these are well understood as coherent states (see e.g. ref. 
|2l|). The coherent states have long been regarded as a close quantum approximation to 
idealized classical states and therefore constitute a natural choice for the preferred basis 
states in a quantum state reduction model. 

The paper is organized as follows. In section [Til we demonstrate the state reduction 
mechanism for the simple case of a quantum harmonic oscillator. By analyzing the quantum 
variance processes we are able to demonstrate that state reduction occurs, and to estimate 
the associated reduction timescale. We also examine how the expectation of energy evolves 
and demonstrate that initial quantum probabilities match with the probabilities of stochastic 
outcomes in a simple example. We conclude the section with some numerical results which 
confirm our analysis. 

In section II I II we extend the formalism to a relativisti c q uantum scalar field. We adopt 



the interaction picture of Tomonaga and Schwinger [22|, |23] to describe a state defined on 
some space-like hyper-surface evolving in a time-like manner. Once we have examined this 
picture in detail, we proceed to demonstrate the reductive properties. We show how this 
mechanism of state reduction for a bosonic field could induce a state reduction to some 
charge state basis in a fermionic field. We end in section [TV] with some concluding remarks. 



II. QUANTUM MECHANICAL HARMONIC OSCILLATOR 

The device we shall use to represent quantum state reduction will be presented for the case 
of (0 + l)-dimensional scalar field theory, i.e. the quantum mechanical harmonic oscillator. 
The commutation relation between position and momentum is given by [x,p] = i. We define 
creation and annihilation operators in the standard way as follows 



^{x — ipuj ip = —iy/~{a — a 

These operators satisfy the commutation relation [a, at] = 1. The Hamiltonian for the 
harmonic oscillator is given by 

H = \p 2 + \u 2 x 2 =u(a?a + l) = u{N + ±), (3) 

where N = a^a is the particle number operator. Units are chosen such that h — 1 for the 
sake of simplicity. 

The Schrodinger equation expressed in differential form is d\ip t ) — —iH\i/) t )dt. We extend 
this in the following way 

d|^> = { [-iH - |AV - a t )a + \\ 2 {a - a t )a t \ dt + X(a - a t )dB t ) |^>, (4) 

where 

at = KMa + a^Nk), ( 5 ) 



4 



and A is a constant parameter of dimension [time]~ l l 2 . Denoting unconditional expecta- 
tion with respect to the physical probability measure P by E [•], the differential dB t is an 
increment of real P-Brownian motion with the properties that E p [di? t ] = 0, idB t ) 2 = dt, 
and increments at different times are independent. Equation can be derived (see 
by first assuming a state evolution equation of the form d\4> t ) = (Cdt + \adX t ) \<f> t ) where 
l^t) — |0t)(0t|0i) _1,/2 and where {X t } is a Q-Brownian motion. The physical measure P 
is related to Q through F(A) = E p [l^] = E^[(^|0*)1a] f° r some event A measurable at 
time t, where 1a — 1 if A is true and otherwise. This choice of physical probability mea- 
sure is the counterpart to the postulate of standard quantum mechanics on the outcomes of 
measurement processes 0. 

Note that since the state evolves according to equation (pO) by the action of only the 
number operator and the annihilation operator, a final state with higher energy than any of 
those states contributing to the initial superposition \ip ) cannot occur. This ensures that 
as long as the initial state has finite energy, subsequent evolved states must also have finite 
energy. 

We proceed by demonstrating that equation (j3J) preserves the norm of a state. Denoting 
\dipt) — dlV't) we have 

d((^|Vt» = (<Wt> + (i>t\dip t ) + <dVt|dVt> 

= (Vtl [M - \\ 2 a\a - a*) + \\ 2 {a ] - a t )a t ] \^ t )dt + (^|A(a f - a t )\^ t )dB t 

\-\H- \\ 2 {a) -ck)a+ \\ 2 {a-a t )a t ] \ijj t )dt+ (ip t \ \(a - a t ) \ip t )&B t 

+ (il) t \\ 2 (a ] - at) (a - a t )\tp t )dt 

= 0. (6) 

For convenience we take the norm of the initial state \ipo) to be unity. Further, we make 
the following definitions for the conditional expectation and conditional variance of some 
operator O with respect to the state \ipt) at time t 

O t =^ t \0\^ t ) and V t ° = {i; t \(0^-0*)(0-OtM), 
and the conditional covariance of two operators O and O' 

v t °>°' = (ip t \(ot-o*t)(o'-o'tM). 

In addition, we define the operator AO t = O — O t . 

Let us first consider the energy of the oscillator. It is straightforward to demonstrate 
that the energy process H t = {ip t \H\ip t ) satisfies the evolution equation 

dH t = -\ 2 ujN t dt + \uj(i) t \(a ] a ] a + a ] aa - 2a Jf aa t )\ip t )dB t . (7) 

By integrating and taking the unconditional expectation we infer that 

= H - X 2 u [ duE p [N u ] . (8) 
Jo 

The second term on the right side is negative semi-definite. Therefore, energy is lost from 
the harmonic oscillator on average at a rate determined by A 2 . We demand that energy 
loss on a macroscopic scale is negligible in order to conform with the energy conservation 



E r [H t ] =H - A 2 E P 



duuN„ 



principle. Taking the typical particle number in the state \ifj t ) to be of order N , we therefore 
require that X 2 uN Q At -C Hq for typical timescales At. Equivalently we may say that A must 
be very small in standard macroscopic units of time. In this limit we have that E p [if t ] ~ H Q , 
or that the expected energy is approximately conserved. In addition, having very small A 
means that for a small number of particles, equation (Tj0) can be accurately approximated by 
Schrodinger's equation. 



A. State reduction 



In order to see how the state reduction mechanism works we consider the stochastic 
processes a t and V t a for the conditional expectation of the annihilation operator and the 
associated conditional variance: 



da t = —iuja t dt — ^\ 2 a t dt + \(ipt\ [(a + a')a — 2a t a\ \^ t )dB t , 
dV t a = -A 2 {(^||Aa,| 2 |^> + |(^| [(0 + 0^)0-20*0] |^)| 2 }dt 
+A(V>t| [(a f - a t )\Aa t \ 2 + \Aa t \ 2 {a - a t )] \^ t )dB t . 

Integrating and taking the unconditional expectation of equation (fTUj) we have 



(9) 
(10) 



E r [V t a ] = V a - A 2 E P 









f duV u 


-A 2 E P 


[I 


Jo 







du\Vj 



a+a'),a 1 2 



Vq - A 2 / dwE p [V£] — A 2 / duE p \Vj 



(a+at),a|2 



(11) 



Since the last two terms on the right side are positive semi-definite, the unconditional ex- 
pectation of the variance of a cannot increase (i.e. V t a is a super-martingale). If we suppose 
that these terms are nonzero then E p [VJ a ] — > for large times and therefore V t a — > i.e. 
the state approaches an a-eigenstate. Otherwise, if for some time t we have E p [V t a ] = 



and E 1 



IK 



(a+a*),a 1 2 



0, then \ip t ) at that time must be an a-eigenstate. Note that the 

second of these two conditions is also satisfied when \ij) t ) is a position eigenstate at time t. 
Since these are composed of an infinite number of infinitesimal energy mode contributions, 
we exclude this possibility. 

In order to estimate the characteristic timescale for state reduction we approximate equa- 
tion (ITTj) by freezing the stochastic terms on the right side at t = 0. In this approximation 
we find 



E p [^ a ] 



ya 
v 



ya 
v 



-A 2 1 + 



\Vr 



(a+at),a 1 2 ' 



ya 
v 



t. 



(12) 



Taking V a ~ V c 



(a+a^ ),a 



O(N ) (corresponding, for example, to a superposition between a 



large excited state and the vacuum state), the reduction timescale for the variance-decreasing 
process can be estimated as 
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This must be small in standard units for macroscopic objects such that macroscopic su- 
perpositions are suppressed. For example, for an oscillator with frequency of order 10 14 s _1 
(corresponding to visible light), if we take N Q = 10 23 and h = 10 _34 Js, then choosing 
A = lO^s- 1 / 2 would lead to energy loss at a rate of 10 13 Js 1 and state reduction on a 
timescale of order 10~ 7 s. For one particle (Nq = 1) energy loss is of order 10~ 36 Js -1 and the 
reduction timescale is 10 16 s (10 9 yrs). 

Once the system enters an a-eigenstate, equation (J2J reduces to 

da t = (-iu; - \\ 2 ) a t dt, (14) 

with solution a t = aoexp{— iut — |A 2 t}. The solution decays on timescale A~ 2 which as 
stated earlier must be very large in standard macroscopic units of time. 

So far we have demonstrated that our state evolution equation (BJ describes state reduc- 
tion to a coherent state on timescale r# given in equation f[T51) . We have also shown that 
coherent states themselves will decay to the vacuum state on a very long timescale A -2 . We 
conclude this subsection by confirming that stochastic probabilities match with quantum 
probabilities for the outcome of a simplified measurement. Let us consider the projection 
operator of a particle number eigenstate P n = \n)(n\. The conditional expectation of the 
projection operator P n>t = (i/; t \P n \i/; t ) obeys the evolution equation 

dP n>t = \ 2 [{n + l)P n+1 , t - nP n>t \dt + \(4> t \{a ] P n + P n a - 2a t P n M)dB t , (15) 

where the terms in square brackets on the right side corresponds to the background de- 
cay mechanism occurring on timescale A -2 . These terms together are small when a given 
wavepacket is sufficiently smoothly varying in n. (For example, a wavepacket centered at 
n = n' with a standard deviation in n of 0(\fnJ), typically has P n> t ~ 0{l/\fnJ) and 
[P n +i,t - Pn,t\ ~ 0(l/n'), resulting in [(n + l)P n +i,t ~ nP n>t ] ~ 0(1). These orders of 
magnitude correspond to a minimum uncertainty coherent state wavepacket.) 

Consider now an initial superposition state |^ ) consisting of the vacuum state |0) and 
some excited coherent state \a ). Suppose further that (0\a t ) ~ 0. We may think of this 
situation as corresponding to a superposition of null and positive readings on some measuring 
device. 

After some time t where tr <t A~ 2 reduction has occurred to a coherent state. This 
may be either the vacuum state or \a t ). The initial quantum probability for registering the 
system in the vacuum state is P va c,o = (^ol-F'vaclV'o) where P vac = |0)(0|. From equation fTl5|) 
we have (upon ignoring the terms in square brackets) 

dP vac , t ~ A(^| (atp vac + P vac a - 2a t P vac ) |^)dP t . (16) 

Now taking the unconditional expectation we have 

P vac ,o ^ E p [P vac , t ] ~ E p [l|^ )= |o>] • (17) 

The final approximation results from the fact that the state at time t is either the vacuum 
state or the approximately orthogonal excited coherent state \ot-t)- This relation tells us 
that the initial standard quantum estimate for the probability of finding the system in the 
vacuum state is equal to the stochastic probability of that outcome occurring in this model. 
The quantum and stochastic probabilities for the other outcome must also be equal. 
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FIG. 1: Conditional expectation of energy. The plot shows five realized paths for an initial state 
corresponding to an equal superposition of two coherent states with expected energies 0.5 and 64.5 
respectively. In the cases where the state reduces to the excited coherent state we note a slow decay 
in energy. This is expected to occur on a timescale of order A -2 ~ 4 in this example ( A = 0.5 and 
to = 1). 



B. Numerical simulations 

In order to confirm the reductive properties we ran a numerical simulation of the quantum 
state evolution. We considered an initial state corresponding to an equal superposition of 
two a-eigenstates with eigenvalues and 8 respectively. We have set the parameters to 
A = 0.5 and uo = 1. This choice means we observe state reduction for small numbers of 
particles with only a small degree of energy loss. Since iVo ~ 32 we estimate the reduction 
timescale by equation (fT3|) to be tr ~ 0.125. The decay timescale is given by A -2 ~ 4. These 
order-of-magnitude estimates are confirmed by figures 1 and 2 which show sample paths for 
the conditional expectation of energy and for the conditional variance in a respectively. 
We see that the state evolves into either one of the two possible coherent states. One of 
these states is the vacuum state, the other corresponds to the (slowly decaying) non- vacuum 
coherent state. 

In addition we have estimated the probabilities of the two possible outcomes by running 
100 sample paths. We find probabilities of 0.47 for the vacuum state and 0.53 for the 
non- vacuum state (the standard deviation of this estimate is 0.1). 



III. RELATIVISTIC QUANTUM FIELD THEORY 



Here we generalize the analysis of the previous section to the case of relativistic quan- 
tum field theory. (For a discussion of the conceptual issues surrounding the formulation of 



relativistic state reduction models, see [14j, [24j, |25|, |26|.) Given that experimental evidence 



conforms to the principle of relativistic invariance it is natural to require this condition of 
our model. This has been a longstanding problem in the field of dynamical state reduction 
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FIG. 2: Conditional variance of the annihilation operator. The sample paths correspond to those 
in figure 1 ( A = 0.5 and u> = 1). 



models. The reason is that while state reduction can be modelled easily enough, by coupling 
a stochastic process to a quantum field we generate energy at an infinite rate. We will resolve 
this issue by coupling only the annihilation operators of the quantum field to the stochastic 
process (as in the case of the harmonic oscillator discussed in the previous section). This 
will ensure that energy cannot be created from the vacuum. 

A natural formulation of relativistic quantum field theory for the consideration of an 



evolving state is the one due to Tomonaga and Schwinger [13j, |22|, |23| . We write the Hamil- 
tonian density at space-time point x in the form H(x) = H Q (x) + H int (x), where H is the 
free field Hamiltonian and H int is an interaction term. Then evolution of the quantum state 
is described by the Tomonaga equation: 



5a(x) 



\*{a)) = H int {x)\*{a)). 



(18) 



The state is defined on some space-like three-surface a, and functional differentiation is 
defined with respect to some point x lying on a. Given two space-like surfaces a and a' 
differing only by some infinitesimal spacetime volume du x at point x (see figure 3) the 
functional derivative can be expressed as 



jljKg)) 

Sa(x) 



lim 



|tt(qQ) - |tt(<7» 
dc^T 



(19) 



Equation (jl8p describes the evolution of the quantum state in terms of incremental time- 
like advancements of individual points on a space-like surface. The operator H int must 
be a scalar quantity in order that equation (I18p has a relativistically invariant form. In 
addition, the constraint [H int (x), H int (x')} = for space-like separated x and x' is imposed 
so that the ordering of points undergoing time-like advancement is irrelevant. We will 
consider the possibility of a definite ordering of all space-time points, allowing us to break 
the commutation constraint. 
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FIG. 3: Evolution between space-like hyper-surfaces a and a'. 



In differential form the Tomonaga equation can be represented as follows 

d x \*(a)) = -iH int (x)\V(a))duj x . (20) 
We proceed by generalizing this equation to a diffusion equation. 

A. Field state diffusion equation 

Previous approaches to modifying Schrodinger field dynamics have generally involved the 
inclusion of a white- noise field term in the Tomonaga equation (see e.g. 0]). Here we opt 
to formulate our model in terms of a Gaussian process. We begin by defining dW x to be an 
increment of some real Q-Brownian motion with mean zero and covariance given by 

E Q [dW4dW4/] = 8 XtX ,dw x . (21) 

We may think of the Gaussian random variable W(a) defined on some surface a and of 
dW x as the incremental difference in W between two surfaces differing by some infinitesimal 
space-time volume at point x. 

We extend the differential Tomonaga equation to include a stochastic term as follows: 

d x \<&(a)) = {-\\ 2 A\x)A{x)du x + \A{x)dW x )\§{a)). (22) 

Here A(x) is a scalar operator to be specified. 

When using the Tomonaga picture, in order to set the initial conditions we must specify 
an initial state on a definite initial space-like surface ov If we then wish to calculate the 
expected state at a later localized region in spacetime, we must specify a final space-like 
surface Oj which includes that region. To describe evolution from the initial state on <7j 
to the final state on Of we could choose any causally ordered set of intermediate space-like 
surfaces (we write a' > a if a' is nowhere in the past of any point on a). Each surface 
will differ by only an incremental spacetime volume du x from its neighboring surfaces in 
the ordering. If evolution of the state from the initial to the final surface is independent of 
the ordering of intermediate surfaces we can say that it is independent of any specific local 
frame. This is true of equation ff22|) provided that [A(x), A(a;')] = [A(:r), Af(x')] = for 
space-like separated x and x'. 
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Assuming the usual rules of Ito calculus we find 

cu(<<t>(^)i<^)>) = 2\{A{x)) <T m<j)m<T))dw x , (23) 

where A(x) = \(A(x) + A\x)), (•)„ = (#(<r)| ■ |#(<r)>, and \^{a)) = |$(a))($(a)|$(a))~l 
is the normalized state. The solution to this equation can be formally written as 

= + 2\ H (A(x)) a ($(a)ma))dW x (24) 

J ai 

= ($(a. J )|$( ( T,))exp|2A^ ' \A(x))AW x -2\ 2 jf ' (A(x))^.| . (25) 

We next introduce the physical measure P such that for a random variable X, measurable 
on surface c/, the P-expectation is given by 



E P [X] = E< 



[*(*/)) ' 
<*fa)l*0*)> . 



(26) 



The physical measure P assigns physical probabilities to possible measurable outcomes. We 
have from equation (124p that P(O) = E p [l] = 1 as required of a probability measure. Also, 
as a consistency check, given the tower law of Q-expectation, it can be shown that the tower 
law of P-expectation also holds: 

E P [X] = E P [E P [A>]]. (27) 

Here a is some surface such that af > a > <Xj, and by conditioning on a we mean that all 
dW x to the past of a are known. It therefore makes no difference for the final outcome if we 
condition on some intermediate surface before taking the expectation at a^. The application 
of equation (1261) therefore provides a consistent way of assigning physical probabilities to 
outcomes. This allows us to describe state evolution in terms of the Q-Brownian motion 
before using the P-measure to determine physical probabilities at the end of the calculation. 

We can also express the state evolution equation directly in terms of a P-Brownian motion 
as follows. First we choose a definite sequence of space-like hyper-surfaces {a} (with <jf > 
a > a,j) to support our state evolution. We then define the process B(a) by the solution to 
the following stochastic equation 

dB x = dW x - 2X(A(x)) a du x . (28) 

Here a is different from its succeeding surface only by some incremental space-time volume 
about x. It can be shown that E p [d-B x ] = and K ¥ [dB x dB x /] = 8 XjX rdu x . Therefore dB x is an 
increment of P-Brownian motion. Finally, writing equation (|22p in terms of the normalized 
state |^(cr)) and the P-Brownian motion dB x we find 

d^(a)> = (a(x, a)du x + (3(x, a)dB x ) |*(<r)), (29) 

where (cf. equation (jlj)) 

oc{x,a) = -\\ 2 {A\x) - (A(x)) a )A(x) + \\\A{x) - (A(x)) a )(A(x)) a , (30) 
f3(x,a) = X(A(x)-(A(x)) a ). (31) 
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In the case where = A^(x')] = for space-like separated x and x', this 

evolution equation must be independent of the choice {cr} by construction. On the other 
hand, if we allow for [A(rr), = [A(x), A^(x')] 7^ 0, the state evolution described by 

equation (j2"2"j) is {<r}-dependent and our choice of sequence {a} must be specific if the model 
is to give unambiguous results. Equation ( 1291) retains its relativistically invariant form so 
that all observers will agree on outcomes. 

We will be forced to choose operators A(x) that do not commute at space-like separation 
and therefore we must specify a fixed sequence of evolving space-like hyper-surfaces. This 
might seem a significant compromise, however, it is not clear that a freedom to choose any 
space-like surface is desirable in a model of quantum state reduction. Consider an entangled 
EPR pair where one particle is measured at a space-like separation from a region where 
we wish to consider the state of the other particle. The state of the unmeasured particle 
depends on whether the surface on which it is defined has the measurement event in its past 
or future. Since we are free to choose this surface, the state of the unmeasured particle is 



ambiguous. As is argued in reference [14J this problem only persists for the state reduction 
timescale so it can be ignored for macroscopic objects. However, it is a difficulty if we intend 
for our state to represent the microscopic world unambiguously. 

If nature were to choose the specific sequence {a} this problem could be avoided. We 
would have no freedom to choose the surface upon which the final state is defined. We 
do not suggest a rule for the choice. We only suggest that relativistic invariance could be 
recovered in expectation by assuming that future space-like surfaces are chosen at random 
from a uniform distribution over the space of all future space-like surfaces. Alternatively, we 
might simply be content to allow our model to break relativistic invariance in its description 
of state reduction. For example, the evolving surfaces could correspond to the constant time 
surfaces in the co-moving frame of the Universe or to a local frame defined by the matter 
content of the state. 

Without a rule for choosing the sequence of surfaces we must quantify the effect of making 
different choices. In fact, we can demonstrate that the imposed ordering of space-time points 
has negligible effect in a perturbative calculation scheme involving the coupling parameter 
A. Given some operator O such that d x O = in the Tomonaga picture, we can use to 
determine the dynamical equation satisfied by its conditional expectation as 

d x (0) a = {a\<r)0 + Oa(a) + ft(a)Ofta)) a duj x + (ft(cr)O + OP(*)) dB x , (32) 

(where the x dependence of a and (3 is assumed). Integrating and taking the unconditional 
expectation we find 



E p [{O) oj ] = (0} ai + E ff 



/ (a\a)0 + Oa(a) + ft (a)0/3(a)) a duj x 



(33) 



Since a ~ 0{\ 2 ) and (3 ~ O(X), we can expand E p perturbatively in A to second 

order by freezing the stochastic state at the initial surface cr*, that is, 



E p [(0) af ] ~ (0) ai + E p [ (at(<Ti)0 + Oa(<Ti) + ft (a t ) 0(3 (a,)) ai du; x 



(34) 



In this approximation, even when the v4(x)-operators do not commute at space-like sepa- 
ration, the result only depends on the choice of initial and final surfaces, and not on any 
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ordering of space-time points within the integrated region. The choice of intermediate sur- 
faces will have no effect. We will use equivalent frozen state approximations in subsequent 
sections. The results will be Lorentz invariant in the sense outlined here. 



We end this subsection by commenting on the "Free Will Theorem" [27j which claims to 



show that relativistic dynamical reduction models are incompatible with the experimenter's 



free will to decide which observable to measure. In subsequent responses [28|, [29j it has 
been argued that the resolution of this conflict can be found in nonlocality (see also 30]). 
Certainly equation (}2"9"]) is explicitly nonlocal through its dependence on the quantum state 
over the entire space- like surface a. However, as pointed out by 't Hooft [3l[, for models of 
this type we should reconsider our notion of "free will". For example, given some definite 
quantum state defined on some initial surface a,, and given some realized B(a) for every 
space-like surface a to the future of <jj, then the future quantum state is determined. This 
future quantum state should describe all matter including the experimenter's behavior. If 
we require free will in this framework, it can only result from an inability to determine the 
precise initial state (3l| . 



B. Scalar field theory 

Having established the covariant form of the theory, we now focus on a particular frame 
with space-like surfaces chosen to be the constant time surfaces. We have 

|d*(t)> = I dJtfC*)) = I dx(a(x)dt + P(x)dB t (x.))\y(t)} } (35) 

</ x ./ x 

with E p [d5 4 (x)d-B t /(x')] = <5 3 (x — x.')5t t t>dt. We use the integration subscript to avoid con- 
fusion over which variables are integrated over. In this frame, time-independent operators 
in the Schrodinger picture are related to time-dependent operators in the Tomonaga picture 
by the unitary transformation 0(t) = exp{iHot}0 exp{— iHot}, where Hq is the free field 
Hamiltonian. 

We consider a real scalar field ip defined in the Tomonaga picture by 

/dp 
. {exp (ip ■ x - itUpt)a(p) + exp (-ip ■ x + iu p t)a\p) \ , (36) 



with free Hamiltonian 



H = J dx {| (d t cp(x)) 2 + \Vip{x) ■ Vip(x) + \m 2 ip\x)} 

= J dp^ p {at(p)a( P ) + |5 3 (0)}, (37) 

where u p = ^/p 2 + m 2 , and the creation and annihilation operators satisfy the canonical 
commutation relations [a(p), a'(p')] = ^ 3 (P — p') an d [ a (p))0(p')] = 0. The positive and 
negative frequency components of the field are given by 

/dp 
7= — ex P (ip ' x - k>p*Mp), (38) 

/dp 
— exp (-ip ■ x + i^ p t)a f (p), (39) 
V 2uj p 



if [X 
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where <p = (p + + ip . We define 



a 



1 \2 
2 



A 2 (<P~ - W)t) ^ + |A 2 (<p + - \{<p) t ) \{v) t (40) 
(3 = \ {^-\^) t ). (41) 

Here (-) t = (^{t) \ ■ \^(t)). The constant parameter A in this model has dimension [time]^ 1 . 
We ignore for now any other possible Hamiltonian interaction terms. Since (p + and do 
not commute, our choice of constant-time surfaces must be considered special. Although, as 
we have seen in the previous section, by using the frozen state approximation, our results 
will be independent of any specific local frame. 



In the same manner as ([6]) we can demonstrate that 

d(*(t)|*(t)) = 0, (42) 

so without loss of generality we may set (^/(t)|\E'(t)) = 1 with the state remaining normalized 
for all time. 

Given some generic operator 0(t) in the Tomonaga picture, we may ask how its condi- 
tional expectation evolves. We find (cf. [l5| ) 



d(0) t = (dO) t + ! dx(a t O + Oa + ftOp) t dt + [ dx^O + O0) t dB t (x) 



(43) 



where dependencies on spatial coordinates are understood. The first term on the right 
side results from the standard unitary evolution of the operator O described by the free 
Hamiltonian. 

Similarly we can write an evolution equation for the conditional variance of an operator. 
Recalling that AOt = O — (O)t and that the conditional variance is given by Vt[0] = 
(|AO t | 2 ) t , we find (again cf. 0) 

dV t [0] = (dO j AO t + AOldO) t + I dx(a t |AO t | 2 + \ AO t \ 2 a + ft\AO t \ 2 (3) t dt 

J X 

- / dx(pW + oip) t (pO + op)tdt 

J X 

+ f dx(/?t|AO t | 2 + \AO t \ 2 p) t dB t (x). (44) 

J X 

Note that the third term on the right side of equation (1441) is negative semi-definite. This 
term is responsible for the variance reduction which we can use to demonstrate state reduc- 
tion (see next subsection). 

We may apply equation f|43|) to the total energy of the quantum field. Ignoring the 
vacuum energy and interactions, this is given by 

H = J dp^ p at(p)a(p). (45) 

We find after some calculation that 

d(H) t = -|A 2 (iV> t dt + / dx(/3t# + HP) t dB t (x), (46) 

J X 
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where 



N 



dpa^(p)a(p) 



(47) 



Integrating and taking the unconditional expectation of the energy process at time t we have 



E v [(H) t ] = (H) - \\ 2 fduE^N)^ 

Jo 



(48) 



(cf. equation (jSJ)). Since (N) t is nonnegative, it follows from ()48l) that energy is lost on 
average as a result of coupling the quantum field to a classical stochastic process. However, 
the energy loss is finite and can be made negligible by an appropriate choice of A. This is to 
be contrasted with some of the previous attempts to construct a relativistic state reduction 
model 13, ll4|, |l5j, where the energy density is seen to increase at an infinite rate. The 
reason that we do not see an infinite rate of energy density creation can be traced back to 
the fact that the classical stochastic process is not coupled to the particle creation operator 
and therefore cannot drive particle creation from the vacuum. 

Stochastic movements in the energy process will cease when the quantum state is an 
eigenstate of the operator ip + . When this occurs, the final term on the right side of equation 
( |4"S|) goes to zero. 

We can approximate equation (J3BD to 0(X 2 ) by freezing the stochastic terms on the right 
side at time t — 0. This gives 



E r [(tf> t ]^(tf)o-!A 2 (iV>ot 



(49) 



This result depends on the initial state and on the integrated region of space-time between 
the initial and final space-like hyper-surfaces. However, no ordering of space-time points is 
required. 

C. Quantum field state reduction 



To see the reductive properties we consider the particle annihilation operator a(p). Using 
equation (j4"51) we find 



d(a(p))t = -ia> p (a(p)) t dt 



(a(p)) t dt + A / dx((v? - (<p) t ) a(p)) t dB t (x). (50) 



Similarly using equation (1441) and taking the unconditional expectation we have 



E ¥ [V t [a(p)}} = V [a(p)] 



2u r 



-E 



d«K[a(p)] 



o 



-A 2 E P 



du / dx|((y> - (<p) t )a(p)) t \ 

I! Jx 

2 rt 



v Q [a( P )] - 



x 2 



2u 



duE p [KKp)]] 



p 



—A duE 



dx|((^- ((p)t)a{p))t 



(51) 
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Again we find that the conditional variance for the annihilation operator is a super- 
martingale. The expected variance decreases with time and the quantum state evolves 
towards an eigenstate of the annihilation operator. If we freeze the stochastic terms on the 
right side of equation (I5T1) we find 

E p [V^[a(p)]] ~ V [a(p)} - ^V [a(p)]t - A 2 jf* du J dx\{{<p - (<p)o) a(p))o| 2 . (52) 



Note that, as in equation (j49j) . the right side is independent of the ordering of space-time 
points and therefore independent of the intermediate space-like hyper-surfaces we have cho- 
sen to support our state evolution. We may estimate the timescale for collapse in the same 
manner as equations f|T2l and (IT31) by taking Vo[a(p)] ~ N (p) = ( a Hp) a (p))o an d 

^dxK( v - Wl)a (p)) ( |^/dp'^Eg2(P), (53) 

from which we find 

TR ~ P/dpWo(p')/(2c p ,)- (54) 

As in the harmonic oscillator case, it is the third term on the right side of equation 
( !5T|) that leads to variance reduction for macroscopic energy scales. The reduction time is 
inversely proportional to the total number of excitations in all modes. This will lead to 
rapid reduction for large scale excitations. Each mode tends towards a coherent state. As 
this occurs, we expect that the field tends towards classical behavior. 



D. Fermionic state reduction 

Here we introduce a fermionic field coupled to our proposed scalar field theory in order 
to consider an induced state reduction in the fermionic sector. To see how this works let us 
set A to zero for now and consider an interaction Hamiltonian of the type 

H- mt {t) = j dxj{x)ip{x). (55) 

J X 

Here j is some Hermitian current operator associated with the fermionic matter field. From 
equation (13"6T) we have 

H int (t) = J dxj(x) J {exp (ip ■ x - iu p t)a(p) + exp (-ip ■ x + i^ p t)a f (p)} 

/dp 
-7= { j f (p, t) exp (-i^pt)a(p) + j(p,t)exp(i^ p t)a t (p)} . (56) 

Furthermore, we can formally solve the Tomonaga equation to find 

|tf(t)) = ex p{"i^ d«tf int («) j|#(0)>. (57) 
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Now suppose that the fermionic state undergoes some spatial transfer of charge such that a 
pulse of current occurs. If the fermionic state is a j-eigenstate, we have 

|*(t)> = exp dp(a(p,t)at(p) - a*(p,t)a(p))W(0)), (58) 
where the complex number a is given by 

a(p, t) = — i / du ^-^ \ exp(iu) p u) , (59) 
Jo a/2w p 

and j(p,t) is the current eigenvalue at time t. Using the commutation relations for the 
creation and annihilation operators, and assuming that the initial tp state is unexcited, we 
find 

a(p')|*(*)> = a(p')exp(/ dp(a(p, t)a\p) - a*(p, t)a(p)) Wo)> 



= exp|y dp(a(p,t)a t (p)-a*(p,t)a(p))j(a(p , ) + a(p , ,t))|*(0)) 

= a(p',f)|¥(f)>. (60) 

The final state is a ^-coherent state with eigenvalue a (cf. section 3.4 in |21j). This demon- 
strates that coherent states in if are associated with j-eigenstates in the matter field. Re- 
duction to a (^-coherent state should therefore induce reduction to a j-eigenstate in the 
fermionic sector. 

It is tempting to associate ip with a gauge field such as the photon field or some pro- 
posed graviton field. The current j would then relate to a conserved charge, e.g. electric 
charge or energy-momentum. Such charge densities are a natural description of macroscopic 
observables. 



IV. CONCLUSIONS 

The key advance of this paper has been to develop an alternative model of state reduction 
in relativistic quantum field theory which does not suffer from the infinite rates of energy 
density increase seen in some previous proposals. We have outlined a model requiring 
just one extra parameter in addition to those of standard quantum theories in order to 
simultaneously describe the quantum behavior of individual excitations and the definite 
behavior of macroscopic objects. 

In our approach, by having no coupling between the classical stochastic field and the 
particle creation operator, we ensure that the evolution equation cannot randomly create 
particles from the vacuum. Our model features only a coupling between the stochastic field 
and the particle annihilation operator. This is appealing for two further reasons. First, 
it leads to a reduction to coherent states. As coherent states saturate the bound of the 
Heisenberg uncertainty relation they make a natural choice as a quantum counterpart to an 
idealized classical state. Second, by applying this mechanism to a bosonic field coupled to a 
fermionic field, we can induce state reduction to some charge density basis in the fermionic 
sector. The model requires the specification of a preferred set of space-like hyper-surfaces 
supporting the time-like state evolution. This breaks relativistic invariance. However, our 
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perturbative calculations show no deviation from relativistic invariance to second order in 
A. 

The ideas presented in this paper could be applied to the photon field or to a pro- 
posed graviton field in order to see state reduction to a conserved electric charge or energy- 
momentum basis in the associated matter fields. Since the model predicts an energy loss 
which could be significant in high-density highly accelerating matter environments, there 
may be the possibility of experimental investigation, e.g., by looking at the decay of high 
intensity electromagnetic waves or through the detection of gravitational waves. 
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